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j^ , Abstract 

■*-» . 
^ ' One of most intriguing problems of heterogeneous crystal nucleation in droplets is its 

strong enhancement in the contact mode (when the foreign particle is presumably in some 
kind of contact with the droplet surface) compared to the immersion mode (particle im- 
mersed in the droplet). Many heterogeneous centers have different nucleation thresholds 
when they act in contact or immersion modes, indicating that the mechanisms may be actu- 
ally different for the different modes. Underlying physical reasons for this enhancement have 
remained largely unclear. In this paper we present a model for the thermodynamic enhance- 
ment of heterogeneous crystal nucleation in the contact mode compared to the immersion 
one. To determine if and how the surface of a liquid droplet can thermodynamically stim- 
ulate its heterogeneous crystallization, we examine crystal nucleation in the immersion and 
contact modes by deriving and comparing with each other the reversible works of formation 
of crystal nuclei in these cases. As a numerical illustration, the proposed model is applied to 
the heterogeneous nucleation of Ih crystals on generic macroscopic foreign particles in water 
droplets at T = 253 K. Our results show that the droplet surface does thermodynamically 
■^- ■ favor the contact mode over the immersion one. Surprisingly, our numerical evaluations 

suggest that the line tension contribution to this enhancement from the contact of three 
water phases (vapor- liquid-crystal) may be of the same order of magnitude as or even larger 
than the surface tension contribution. 
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1 Introduction 

The size, composition, and phases of aerosol and cloud particles affect the radiative and chemical 
properties 1,2 of clouds and hence have a great impact on Earth's climate as a whole. On the 
other hand, the composition, size, and phases of atmospheric particles are determined by the 
rate at and mode in which these particles form and evolve. 2-4 

Water constitutes an overwhelmingly dominant chemical species that participates in atmo- 
spheric processes. Consequently, great importance is attributed to studying aqueous aerosols 
and cloud droplets as well as their phase transformations. In a number of important cases atmo- 
spheric particles appear to freeze homogeneously 4-6 For example, the conversion of supercooled 
water droplets into ice at temperatures below about -30°C is known to occur homogeneously, 
mainly because the concentrations of the observed ice particles in the clouds often exceed the 
number densities of preexisting particles capable of nucleating ice. 4,5 Also, it has been suggested 
that aqueous nitric acid-containing cloud droplets in the polar stratosphere freeze into nitric 
acid hydrates via homogeneous nucleation. 6 Understanding how nitric acid clouds form and 
grow in the stratosphere is a topic of current interest because such clouds participate in the 
heterogeneous chemistry that leads to springtime ozone depletion over the polar regions. 2 

However, most phase transformations in aqueous cloud droplets occur as a result of hetero- 
geneous nucleation on preexisting macroscopic particles, macromolecules, or even ions. 3 Hetero- 
geneous nucleation of ice on a microscopic foreign particle can be considered as the adsorption 
of water molecules on a substrate which serves as a template. If the water molecules adsorb in 
a configuration close enough to the crystalline structure of ice, then the energy barrier between 
phases is substantially reduced. Recent work on the heterogeneous nucleation of ice in the 
atmosphere is motivated by the evidence, primarily from modeling studies, that heterogeneous 
freezing may significantly impact the radiative properties, both in the visible and infrared, of 



cirrus clouds. The leading candidates for heterogeneous nucleating centers are the mineral dusts 
(fly ash and metallic particles) and emissions from aircraft, primarily soot. 7-9 Interest in cirrus 
clouds motivated several laboratory studies as well. 10 ' 11 They showed that the presence of var- 
ious foreign inclusions shifts the apparent freezing temperature of droplets upward by as much 
as 10°C. 

Most investigators targeted particulates as the primary heterogeneous ice nucleating centers 
in the atmosphere. Recently, however, increasing attention is payed to the role of films of high- 
molecular-weight organic compounds located on droplets. Such compounds are emitted into 
the atmosphere, especially in regions that are influenced by biomass burning. 12 It was reported, 
for example, that the films of long-chain alcohols and some other organic species can catalyze 
ice nucleation in droplets at a supercooling of only 1°C. 13 ' 14 

So far, the physical mechanism underlying heterogeneous crystal nucleation in droplets 
remains rather obscure. 2 As an additional mystery, many heterogeneous centers have different 
nucleation thresholds when they act in different modes - contact or immersion, indicating that 
the mechanisms may be actually different for the different modes. In the contact mode, the 
ice-nucleating particle contacts water droplet, i.e., touches or intersects its surface, whereas 
in the immersion mode the particle is immersed in the water droplet (Figure l). 2 ' 15 The same 
particle tends to trigger the freezing of an supercooled water droplet at a higher temperature 
in the contact mode than in the immersion one. 2,15,16 

The cause of this enhancement is unknown, but it provides a hint that the water surface could 
be of special interest in ice nucleation. Several investigators have put forward conjectures on 
the mechanism of contact nucleation, all of which depend on the contact of a particle impinging 
upon the droplet surface from air. 2 One of the hypothesis is based on the partial solubility 
of small solid particles whereby active sites at the surface of a particle are subject to erosion 



after it becomes immersed in water. 17 ' 18 Another hypothesis suggests 19 that only those particles 
enhance nucleation in the contact mode which exhibit a strong affinity for water. During the 
initial contact with the droplet (before the equilibrium adsorption is achieved) such particles 
might strongly lower the free energy barrier to ice nucleation at its surface. Another interesting 
explanation 20 suggests that the contact mode enhancement of crystal nucleation is due to the 
mechanically forced rapid spreading of water along the hydrophobic solid surface which forces 
its local wetting and thereby temporarily creates local high interface-energy zones increasing 
the probability of crystal nucleation. While acceptable for some particular cases, all those 
explanations have some inconsistencies and limitations, and so far no rigorous (and general 
enough) theoretical model of this phenomenon has been proposed. 

As a related problem, recently a thermodynamic theory was developed 21 ' 22 that prescribes 
the condition under which the surface of a droplet can stimulate homogeneous crystal nucleation 
therein so that the homogeneous formation of a crystal nucleus with one of its facets at the 
droplet surface (surface-stimulated mode) is thermodynamically favored over its formation with 
all the facets within the liquid phase (volume-based mode). For both unary and multicomponent 
droplets the inequality coincides with the condition for the partial wettability of at least one of 
the facets of a crystal nucleus by its own melt. 23 This effect was experimentally observed for 
several systems, 24 ' 25 including water-ice 26 at temperatures at or below 0°C. 

Clearly, the mode of crystal nucleation is most likely determined by both thermodynamic 
and kinetic factors. However, the partial wettability of a solid by its melt may help to explain 
why, in molecular dynamics simulations of various kinds of supercooled liquid droplets 2 ,28 the 
crystal nuclei appear preferentially close to the surface. Since smaller droplets have a higher 
surface-to-volume ratio, the per-droplet nucleation rates in small droplets tend to be higher 
than in the bulk. Hence it is experimentally easier to observe the crystallization of aerosols, 



having a large collective surface area, than those having a large volume. Recent experiments 16 
on the heterogeneous freezing of water droplets in both immersion and contact modes have 
also provided evidence that the rate of crystal nucleation in the contact mode is much higher 
because the droplet surface may stimulate heterogeneous crystal nucleation in a way similar to 
the enhancement of the homogeneous process. 

In this paper we extend the approach, previously developed in refs.21,22, to heterogeneous 
crystal nucleation on a solid particle (in both immersion and contact modes) and present a 
thermodynamic model thereof. Our thermodynamic analysis suggests that, indeed, the droplet 
surface can thermodynamically enhance crystal nucleation in the contact mode compared to 
the immersion mode. Whether this occurs or not for a particular foreign particle is determined, 
however, by the interplay between five surface tensions and four line tensions involved in this 
process. 

The paper is structured as follows. In section 2 we derive and compare with each other the 
expressions for the free energy of heterogeneous formation of a crystal nucleus on a solid (say, 
dust) particle in the immersion and contact modes. For the sake of simplicity, in this work we 
consider only unary systems, i.e., pure water droplets, but the generalization to multicomponent 
droplets can be carried out as well. Only one kind of foreign nucleating centers is considered, 
namely, those completely wettable by water. Numerical predictions and possible experimental 
verification of the model are discussed in Section 3. The results and conclusions are summarized 
in section 4. 



2 Free energy of heterogeneous formation of crystal nuclei in 
contact and immersion modes 

To determine if and how the surface of a liquid droplet can thermodynamically stimulate its 
heterogeneous crystallization, it is necessary to consider the formation of a crystal cluster in 
the two modes (Figure 1). In the "immersion" mode, the crystal cluster is formed with one of 
its facets on a foreign particle that is completely immersed in a liquid droplet; all other crystal 
facets interface the liquid. In the "contact" mode, the foreign particle touches (i.e., is in contact 
with) the droplet surface and the cluster forms with one of the crystal facets on the particle 
(as in the immersion mode), another facet at the liquid- vapor interface, and all other facets 
making the "crystal-liquid" interface. In these two cases the reversible works of formation of 
a crystal nucleus (critical cluster) should be derived and compared with each other. This can 
be carried out in the framework of the classical nucleation theory (CNT) for both unary and 
multicomponent droplets. In this paper we consider the crystallization of unary droplets. 

The droplet surface can incur some deformation if its crystallization is initiated at its surface. 
The thermodynamic analysis of the case where the crystallization begins at a droplet surface 
can be considerably more complicated when compared to the case where it forms at the surface 
of a bulk liquid. However, one can show 21 that if a 1 /7rR 2 <C 1 (where a 1 is the surface area of 
the crystal facet interfacing the vapor and R is the droplet radius), the formation of a crystal 
at a droplet surface can be considered as crystallization at the surface of a bulk liquid. Under 
conditions relevant to the freezing of atmospheric droplets, crystal nuclei are usually of sub- or 
nanometer size, while the droplets themselves are in submicron to micrometer size range, i.e., 
the above condition is well satisfied. Since the analysis of the freezing of atmospheric droplets is 
our ultimate goal, one can assume the droplet surface to be flat and thus avoid the complexity 



of taking into account the droplet deformation upon freezing. Besides, heterogeneous particles 
serving as nucleating centers can be considered as macroscopic particle of linear sizes much 
greater than the crystal nuclei hence the part of its surface on which the crystal nucleus forms 
can be considered to be flat as well. 

Let us consider a single-component bulk liquid. A macroscopic heterogeneous particle is 
either completely immersed in the liquid or in contact with the liquid-vapor interface. Crystal- 
lization will take place in this liquid if it is in a metastable (supercooled) state. The reversible 
work of crystal formation, W, can be found as the difference between Xft n , the appropriate 
thermodynamic potential of the system in its final state (liquid+crystal) , and Aj n , the same 
potential in its initial state (liquid): W = X^ n — Xi n . Since the density of the liquid may be 
different from that of the solid, the volume of the liquid may change upon crystallization if the 
process is not constrained to be conducted at constant volume. In such a case, strictly speak- 
ing, one cannot calculate W as the difference in the Helmholtz free energies since the volume 
work that the entire system exchanges with the environment should not be regarded as work 
involved in the formation of a local nucleus. As an approximation, the use of the Helmholtz free 
energy is still acceptable since, in the thermodynamic limit, the change in the total volume of 
the system is usually negligible. A better choice for the thermodynamic potential is the Gibbs 
free energy if the system is in contact with a pressure reservoir (since the unconnected volume 
work exchanged with the environment is automatically removed from the Gibbs free energy). 
However, 29 in the thermodynamic limit, the use of either the Gibbs or Helmholtz free energy 
or grand thermodynamic potential is acceptable for the evaluation of W . 

Neglecting the density difference between liquid and solid phases and assuming the crystal- 
lization process to be isothermal, one can say that the total volume, the temperature, and the 
number of molecules in the system, respectively, will be constant. Thus the reversible work of 



formation of a crystal embryo can be evaluated as the difference between Ffi n , the Helmholtz 
free energy of the system in its final state (liquid+crystal+foreign particle), and Fi n , in its 
initial state (liquid+foreign particle): 

W = F fin - F in . (1) 

2.1 Foreign particle completely immersed in the liquid 

Consider a bulk liquid in a container whose upper surface is in contact with the vapor phase 
of constant pressure and temperature. A macroscopic foreign particle is completely immersed 
in this liquid. Clearly, for this system to be in mechanical and thermodynamic equilibrium, 
the particle must be completely wettable by the liquid. Upon sufficient supercooling, a crystal 
nucleus may form heterogeneously with one of its facets on the foreign particle. The crystal is 
considered to be of arbitrary shape with A facets (Figure 2). We will assign the subscript "A" 
to the facet which is in contact with the foreign particle. (Figure 3) 

Let us introduce the superscripts a, f3, 7 and 5 to denote quantities in the liquid, vapor, 
crystal nucleus, and foreign particle, respectively. Double superscripts will denote quantities at 
the corresponding interfaces, and triple superscripts at the corresponding three-phase contact 
lines. The surface area and surface tension of facet i (i = 1, ..., A) will be denoted by Ai and <7j, 
respectively. (Anisotropic interfacial free energies are believed to be particularly important in 
determining the character of the nucleation process.) Hereafter, we adopt the definition of the 
surface tension of a solid, a s , as given in chapter 17 of ref.19. Namely, a = f + ^\ Tj/x^, 
where /', T, and fi' are the surface free energy per unit area, adsorption, and chemical potential 
of component i, all attributed to the dividing surface between solid and fluid. In the following, 
we will neglect the adsorption at the solid-fluid interfaces. Thus, by definition, the surface 
tension of the solid will be equal to the surface free energy per unit area. 



Let us denote the number of molecules in the crystal cluster by v. Neglecting the density 
change upon freezing and assuming the equality of pressures in the vapor and liquid, P a = P", 
the reversible work of heterogeneous formation of the crystal (with its facet A on the foreign 
particle) is given by the expression 

A-l 
W imm = zy [ /U 7(p7 ; T ) - fi a (P a , T)] - V"< {P"< -P a ) + ^ fff 7 ^ 7 + uf A\ & - <j a& ' A\ & + T^ & U*^ , 

j=l 

(2) 
where n, P, V, and T are the chemical potential, pressure, volume, and temperature, respec- 
tively, and r is the line tension associated with a three-phase contact line 30 of length L. 

The necessary and sufficient conditions for the equilibrium shape (known as the Wulff form) 
of the crystal are represented by a series of equalities referred to as Wulff 's relations (see, e.g., 
ref.23), which can be regarded as a series of equilibrium conditions on the crystal "edges" formed 
by adjacent facets. For example, on the edge between homogeneously formed facets i and i + 1 
the equilibrium condition is 

ay ay 

a i- = ^ (i = l,...,A), (3) 

where hi is the distance from facet i to a point O within the crystal (see Figure 2) resulting 
from the Wulff construction. 23 

In the case when one of the facets (facet A) is the crystal-vapor interface while all the others 
lie within the liquid phase (see Figure 3), the equilibrium conditions on the edges formed by 
this facet with the adjacent ones (hereafter marked by a subscript j) are given by 

"T <-?-<* (4) 

hj i, x ' v 

Note that the height of the A-th pyramid (constructed with the base on facet A and with the 
apex at point O of the Wulff crystal) will differ from that with all of the facets in the liquid. 
Thus, the shape of the crystal will differ from that in which all facets are in contact with the 



liquid. For this case, Wulff's relations take the form 



aZ 1 _ aZ 1 _ _ oT ~ °°* 



(5) 



hi h 2 h\ 

(see also refs.21,22). In the above consideration, it is assumed that the mechanical effects within 
the crystal (e.g., stresses) reduce to an isotropic pressure P 1 . In this case 21 

2a ai 2(rr' yS - rr aS ) 

pl_ pa = J^ (i = l,...A-l), p7_p« = AA^ L, ( 6 ) 

hi h\ 

Equation (6) applied to the crystal is the equivalent of Laplace's equation applied to liquid. 
Thus, just as for a droplet, one can expect to find a high pressure within a small crystal. It is 
this pressure that is the cause of the increase in the chemical potential within the crystal. 

Note that the line tension contributions to the free energy of crystal formation were omitted 
in the model for homogeneous crystal nucleation in the surface-stimulated mode 21 ' 22 because 
they were assumed to be negligible compared to the volume and surface contributions. However, 
this assumption may no longer be valid for heterogeneous crystal nucleation because the nucleus 
is now much smaller (compared to the homogeneously formed one) and hence the contributions 
of three-phase contact lines can be more important. 31 ~ 33 

Making use of equations (5) and (6), one can rewrite eq.(2) as 

A-l 
W imm = (/ | |i 7(pa ) T ) _ ^(p«, T)] + £ (jf 7 ^ 7 + uf Af - O a& A\ 5 + T^ S L a ^ S . (7) 

i=l 

In this equation, the first term represents the excess Gibbs free energy of the molecules in the 
crystal compared to their Gibbs free energy in the liquid state. This term is related to the 
enthalpy of fusion Ah by (see, e.g., ref. 23) 

f T dT' 

^(P a ,T)-^(P a ,T) = - Ah—, (8) 

J T 1 

where Tq is the melting temperature of the bulk solid (T < To), and Ah < 0. 
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If the supercooling T — Tq is not too large or, alternatively, if in the temperature range 
between T and To the enthalpy of fusion does not change significantly, eq.(8) takes the form 

^(P a ,T)-fi a (P a ,T) = -A/ilne. (9) 

with B = T/Tq. Thus, one can rewrite eq.(7) in the following form 

A-l 

/A/i In 9 + J2 °V A V + °t A f - ° aSA f + r a ^L a ^ s . (10) 



j=l 



By definition, the critical crystal (i.e., nucleus) is in equilibrium with the surrounding melt. 
For such a crystal the first term in eq.(2) vanishes. On the other hand, for a crystal with one 
of its facets being a crystal- foreign particle interface, and the others interfaced with the liquid, 
one can show that 

yi^pi _ P «) = 2 fj^ afAri + afAf - a a5 Al s ) . (11) 

(This equality can be derived by representing V 1 as the sum g J2i=i hiAi of the volumes of 
A pyramids with their bases at the crystal facets and their apexes at point O. The difference 
pi _ p a f or every term in this sum is replaced by the RHS of the corresponding equality in 
eq.(6).) Substituting eq.(ll) into eq.(2), one can thus obtain the following expression for the 
reversible work W* mm of formation of a critical crystal: 

W imm = -(j2 °V A V + °-f A f ~ ° a5A l 5 ) + T^L^ 5 , (12) 



or, alternatively, 



yyimm _ _y7(p7 _ p°<\ _|_ T ay8 pa-yS Qg\ 



(hereafter the subscript "*" indicates the quantities for the nucleus; it is omitted on the RHS 
of expressions for Wl mm to avoid the overcrowding of indices) 

Clearly, in the atmosphere the crystal cluster forms not in the bulk liquid, but within a 
liquid droplet (see Figure 1) which is itself surrounded by a vapor phase. The reasoning here is 
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almost identical to the preceding if, again, we neglect the density difference between the liquid 
and crystal phases. One can easily show that all above equations, starting with eq.(2) and 
including eqs.(12) and (13) for the reversible work Wl mm of formation of the critical crystal, 
remain valid except that the pressures in the liquid and vapor phases are not equal but are 
related by the Laplace equation P a = P™ + 2a a ^ /R, with R being the radius of the droplet 
(assumed to remain constant during freezing). 

2.2 Foreign particle in contact with the liquid-vapor interface 

Now let us consider a foreign particle which is not immersed in a bulk liquid but is in some 
kind of contact with the liquid- vapor interface (Figures 1 and 4), with the vapor phase being 
at constant pressure and temperature. (Clearly, the particle must be completely wettable by 
the liquid in order for the same particle to be able to be in mechanical and thermodynamic 
equilibrium in both immersion and contact modes.) Now, upon sufficient supercooling, a crystal 
nucleus may form heterogeneously with one of its facets (marked with the subscript "A" ) on the 
foreign particle and another one at the vapor-liquid interface. The latter facet will be marked 
with the subscript A — 1. All the other A — 2 facets lie within the liquid phase. 

Again, neglecting the density change upon freezing and the equality of pressures in the 
vapor and liquid, P a = P@, the reversible work of heterogeneous formation of the crystal with 
its facet A on the foreign particle, the facet A — 1 interfacing vapor (i.e., in the contact mode) 
and all the others within the liquid phase will be given by the expression 

A-2 



W con = u{^{P\T) - fi a (P a ,T)} - V'^P"' -P a ) + J2 °V A V + ff A-l^A-l " ^^a-i + 



Considering the contact mode, the prime will indicate quantities whereof the values may differ 
from those in the immersion mode. The equilibrium shape of the crystal (the Wulff form) is 
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again determined by a series of equilibrium conditions on the crystal "edges" formed by adjacent 
facets. For example, the equilibrium conditions on the edges formed by the facets A — 1 and A 
with the adjacent ones are given by 

(T al (T® 1 -CT a P rr a ~< rr 1& rr aS 

t^^tr- 0#A) ' t=~^ < ^ A - 1) ' (I5) 

where j and k mark the facets adjacent to facets A — 1 and A, respectively, and primes will 
hereafter mark quantities for facet A — 1 at the droplet surface. Note again that the height of 
the A — 1-th pyramid (constructed with the base on facet A — 1 and with the apex at point O 
of the Wulff crystal) will differ from that with all of the facets (except for facet A, see Figure 
3) in the liquid. Thus, the shape of the crystal will differ from that formed heterogeneously in 
the immersion mode (i.e., when all facets, except for facet A, are in contact with the liquid. For 
this case, Wulff' s relations take the form 



(16) 



°i _ °2 __ _Vx-i~ a aX'-a 



hi h 2 h' x _ x h 

Consequently, eq.(6) (the equivalent of Laplace's equation applied to crystals) becomes 

2a ai 2(a^,-a al3 ) 2(af - a aS 



hi hx-i h\ 

(17) 

Making use of equations (16) and (17), one can rewrite eq.(14) as 

A-2 

w con = v[^(p a ,T) - fi a (p a ,T)] + J2 °V A V + °l lA l\ - ^ A xii + °f A f - ° a&A f + 



Furthermore, using eq.(9) one can represent eq.(18) in the following form 

A-2 

W con = -uAh In 6 + Y, °V A V + ^A-i A \-x ~ ^ A x-i + ef A f ~ ^^x + 

i=l 
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For a crystal with one of its facets being a solid-vapor interface, and the others interfaced 
with the liquid, one can show that 

y'7( P 7 _ p«) = | j £ af 7 A° 7 + afZ^f - a^A^ + af^f - a a<5 ^ j , (20) 

which makes it possible to represent the reversible work Wl of formation of a critical crystal by 
the expression 

Wr = \ (E °V^V + ^a Va-i - ^*X-i + <>?*? - * aS Af) + 

T al3~ f j^al3~f _|_ ^/^ _ r a/3<5 \ £/?7<5 _|_ T a~ f & -^ a~ f 8 ^l) 

or, alternatively, as 

Equations (21) and (22) are similar to eqs.(12) and (13) which apply to heterogeneous crystal 
nucleation in the immersion mode. Along with eq.(13), equation (22) will be most useful in 
later discussions. 

The reversible works of heterogeneous formation of crystal nuclei in the immersion and 
contact modes can now be compared. The difference between the internal pressure of the nucleus 
and the external pressure does not depend on whether the nucleus forms in the immersed mode 
(let us denote it by (P 7 — p a ) jmm ) or in the contact mode (denoted by (Pj — p Q ) con ). Indeed, 
by using equation (9) and the equilibrium condition for the nucleus, namely 

^(P\T)-fi a {P a ,T)=0, (23) 

while assuming the crystal to be incompressible, one can show that the difference P 7 — P a for 
the nucleus, in both cases, is determined by the supercooling of the liquid, so that 

Ah 

(PJ - P a ) con = (PJ - p a f mm = In 6, (24) 

v 
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where v is the volume per molecule in the crystal phase. The first equality in eq.(24) is equivalent 

to 

Pi _ a/3 "7 

a X-i a °A-i 



h\i h X -i ' 



from which it follows that 



a P~1_ _ a af3 

h'x-i = X ' 1 aj h X -i- (25) 

a X-l 

On the other hand, /ij = h!, L for i = 1, . . . , A — 2, A, by virtue of eqs.(6), (17), and (24). This 
means that the Wulff shape of the crystal, in the contact mode, is obtained by simply changing 
the height of the A — 1-th pyramid of the Wulff crystal in the immersion mode. It is hence clear 
that if a^-a^ <a x n _ 1 , then 

^_i < fr*-i => V? < VJ. (26) 

According to eqs.(13) and (22), 
\Y con — W imm = -(V' 1 — V 1 )(P 1 — P a ) + T a 7 S (l,' a ~' S _ L a ~t s \ _|_ ( T /375 _ T af35\j^'Pi5 _|_ T a/3~/ j^ajS 

(27) 
Because of eqs.(25) and (26), if 

a^-a^<al\ (28) 

then the first term on the RHS of eq.(27) is negative. 

In the case of homogeneous crystal nucleation the line tension contributions to the free 
energy of crystal nucleus formation are either negligible or non-existent. 21 ' 22 It allowed one to 
conclude 21 ' 22 that if the condition in eq.(28) is fulfilled, it is thermodynamically more favorable 
for the crystal nucleus to form with its facet A at the surface rather than within the liquid. In- 
equality (28) coincides with the condition of partial wettability of the A-th facet of the crystal by 
its own liquid phase. 23 This effect has been experimentally observed for water-ice 26 at temper- 
atures at or below 0°C. In those experiments, 26 when air was added to water vapor the partial 
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wetting of ice by water transformed into complete wetting, but only for some orientations. 
Besides, the wettability of solids by fluids usually decreases with decreasing temperature. ' 35 
Since the freezing of atmospheric water drops always occurs at temperatures far below 0°C, one 
can expect the partial wettability of at least some facets of water crystals even in the presence 
of air. Furthermore, according to Cahn, 36 perfect wetting of a solid by a liquid away from the 
critical point is not generally observed, i.e. the condition in eq.(28) should be fulfilled for most 
substances. In Cahn's theory, the general restrictions on the solid phase are that its surface 
is sharp on an atomic scale and interactions between surface and fluid are sufficiently short- 
range. Therefore, that theory can be also applied to the case where the temperature is far 
below the fluid critical point and the solid is of the same chemical nature as the fluid phases. If 
the temperature approaches the fluid critical temperature, Cahn's theory becomes inapplicable. 
However, the temperatures involved in crystallization are usually far below the critical point. 
All these combined with eq.(28) helps to explain why, in molecular dynamics simulation studies, 
crystallization begins at or near a surface, and why it is easier, experimentally, to observe the 
homogeneous crystallization of aerosols than that of the corresponding bulk liquid. 

However, the presence of the line tension contribution on the RHS of eq.(27) makes it 
impossible to draw unambiguous conclusions concerning the difference W£ on — W£ mm for het- 
erogeneous crystal nucleation even when inequality (28) is fulfilled. Although the first term on 
the RHS of eq.(27) is negative and gives rise to the thermodynamic propensity of the crystal 
nucleus to form with the facet A — 1 at the droplet surface, the line tension contributions can be 
both negative and positive because any of the line tensions involved can be either negative or 
positive. 31-33 Moreover, the sign of the line tension may change depending on the temperature 
at which the crystallization occurs. 
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3 Numerical evaluations and discussions 

To illustrate the above theory by numerical evaluations, consider first the homogeneous freezing 
of water droplets (surrounded by water vapor in air) at around Th m = 233 K (i.e., about — 40°C). 
The estimates for homogeneous crystal nucleation will serve as a reference point to obtain some 
estimates of the relative importance of the line tension contributions on the RHS of eq.(27) for 

ytfcon TjTimm 

'" * "' * 

As reported by Defay et a/. 23 , the rate of homogeneous crystal nucleation in bulk supercooled 
water at this temperature is 7x 10 12 cto _3 s _1 , with the nucleation barrier height W* = 45 A;T/j m , 
the average (over all crystal facets) surface tension of liquid-solid (water-ice) interface a" 7 being 
about 20 dyn/cm (Table 18.1 in ref.23). The surface tensions of liquid- vapor and solid- vapor 
(ice- water vapor) interfaces at Thm = 233 K will be taken to be a a/3 = 88 dyn/cm and <r^ 7 = 103 
dyn/cm, respectively. All these values of a" 7 , a a/3 , and u^ 7 are consistent with the data provided 
in ref.2. 

Let us assume that only the basal facets of the hexagonal ice crystal is partially wettable by 
water at Th m = 233 K. The height of the basal pyramid of the crystal cluster will be denoted by 
hb when the basal facet is at the droplet surface (hereafter a tilde will mark quantities for this 
case) and by h\, when the entire crystal is immersed in the droplet. The corresponding works of 
homogeneous formation of crystal nuclei will be denoted as W* and W*. According to eqs.(31) 
and (32) of ref.21 and eqs.(15) and (16) of ref.22, one can obtain for crystal nuclei: /i* = 2h b ~ 
21 x 10~ 8 cm, K = h + h b ~ 18 x 10~ 8 cm, W* ~ 39.3 kT hm , AW* = W*-W*zi -5.7 kT hm 

These evaluations are for homogeneous crystal nucleation in the volume-based vs surface- 
stimulated modes 21 ' 22 which are equivalent to the immersion and contact modes, respectively, 
of the heterogeneous crystal nucleation. At any particular temperature, the critical crystal of 
heterogeneous nucleation is much smaller than for homogeneous nucleation. On the other hand, 
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the size of the nucleus increases with increasing temperature (i.e., decreasing supercooling). 
Thus, one can expect that for any foreign particle there exists a temperature 233 K < T^t < 273 
K such that the linear size of the crystal nucleus (hence the number of molecules therein, v c ) for 
heterogeneous nucleation is comparable to that estimated above for homogeneous nucleation at 
Th m = 233 K. The energy unit ksT (ks is the Boltzmann constant) varies from 3.2 x 10 -14 erg 
to 3.8 x 1(T 14 erg, i.e, by about 10 %. 

Let us assume, that for a selected foreign particle the temperature T^t — 253 K with the 
thermal energy unit ksT^t — 3.5 x 10~ 14 erg. At this temperature, the first term on the RHS 
of eq.(27) (hereafter referred to as the "surface-stimulation term") can be roughly assumed to 
be equal to AW* because both quantities represent the surface contribution to the difference 
between the free energy of nucleus formation in the surface-stimulated and volume-based modes 
(for heterogeneous and homogeneous nucleation, respectively). Thus, according to the above 

estimates, 

\(W-V:)(P1-P«) ^ AW. „ T hm ^ 

kbTm k B T ht T ht 

As mentioned above, this contribution is negative if facet A — 1 (formed at the liquid-vapor 

interface) is only partially wettable by its melt (i.e., water), which is the case with the basal 

facet of the crystals of hexagonal ice. Consequently, the droplet surface always makes the 

contact mode of heterogeneous crystal nucleation in water droplets thermodynamically more 

favorable than the immersion mode, regardless of the nature of the foreign particle. 

It is virtually impossible to provide general unambiguous estimates for the line tension 

contributions to Wf on — Wl mm in eq.(27). Indeed, the line tension is notorious not only for the 

lack of reliable experimental data but (mostly) for its ability of being both negative and positive 

and take values in the range from 10 _1 to 10 -5 erg (see, e.g., refs. 31-33). Nevertheless, some 

estimates can provide useful insight into the problem of "contact mode vs immersion mode" of 
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heterogeneous crystal nucleation. 

In the second term on the RHS of eq.(27), the difference L a i° —L a i° represents the difference 
between the lengths of the "liquid-crystal-foreign particle" contact line in the contact and 
immersion modes. Clearly, this difference is negative. Considering, as above, that it is the 
basal facet of the hexagonal ice crystal which forms at the liquid-vapor interface (with one of 
the six prismal facets formed on the foreign particle) , one can conclude that L ° 7 — L" 7 ~ 
-a = -8.65 x irr 8 cm for q = 2 and L' a ~< 5 - L a ~< s « -a = -13.8 x 1(T 8 cm for q = 0.5. As 
for the line tension r* 7 , its sign can be expected to be positive, 31-33 but we are not aware of 
any experimental or theoretical data reported for "foreign particle-crystal-vapor" three-phase 
contact regions. Assuming that t" 7 can be anywhere in the range from 1CP 1 erg to 10 -5 erg, 
it is still most likely to be closer to 1CP 5 than to 1CP 1 erg because two out of three phases 
in contact are solid phases involving little inhomogeneities of density profiles in the contact 
region. Thus, one can cautiously suggest that: a) this three-phase contact line impedes the 
heterogeneous crystal nucleation in the immersion mode vs contact mode; b) possible values of 
the term T a ~- ' (L Q7 — L° 7 )//cbT may be somewhere in the range from —100 to —10. 

The third term on the RHS of eq.(27) is due to the three-phase contact line "liquid- vapor- 
foreign particle". The length of this line, L^ 7 , is equal to — (L a7 — L° 7 ), evaluated in the 
above paragraph, i.e., L^ 7 ~ 8.65 x 10 -8 cm for q = 2 and 13.8 x 10~ 8 cm for q = 0.5. 
Further, the density inhomogeneities in the "foreign particle-crystal-vapor" contact region can 
be expected to be negligible compared to those in the "foreign particle-crystal-liquid" , "liquid- 
vapor-foreign particle", or "liquid- vapor- liquid" . Therefore, one can consider that the line 
tension r^ 7 is negligible when compared to r a ^ , so that the third term becomes (r^ 7 — 
r a/3 )L^ 7 ~ — T a P°L@' 1 ° . Depending on the wettability of the foreign particle by liquid water 
in the water vapor, r a/3 can be positive as well as negative. It was noted, however, that for the 
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same foreign particle to be able to serve as an equilibrium nucleating center in both immersion 
and contact modes it has to be completely wettable by watter. Thus, one can suggest 27-29 that 
the line tension T a ^ < with its absolute value closer to 10 -5 than lCT 1 . Besides, the third 
term on the RHS of eq.(27) can be expected to provide a contribution to W£ on — Wl mm which 
is close to the contribution from the second term in absolute value and has an opposite sign. 
The approximate compensation of the second and third terms can thus be expected. 

The last term on the RHS of eq.(27) is due to the three-phase contact line "liquid- vapor- 
crystal" . Considering again the basal facet of the hexagonal ice crystal forms at the liquid- vapor 
interface (with one of the six prismal facets on the foreign particle) , the length of this contact 
line is approximately L a ^ « 5a, that is, L a ^ « 40.2 x 10~ 8 cm for q = 2 and L a ^ « 69.0 x 1(T 8 
cm for q = 0.5. As the basal facet of an Ih crystal is partially wettable by liquid water with 
the contact angle (measured inside the liquid phase) less than 7r/2, one can consider r a/37 to be 
negative. 31 ' 32 Even assuming for the value of r a/37 the smallest experimentally reported order 
of magnitude, 10 -5 erg, one can conclude that: a) this contact line significantly enhances the 
contact mode of heterogeneous crystal nucleation compared to the immersion mode; b) the 
absolute value of the line tension contribution to W£ on — W£ mm from the "vapor-liquid-ice" 
contact line is at least by one order of magnitude greater than that of the surface-stimulation 
term (first term on the RHS of eq.(27)). Thus, this line tension contribution to W£ on — W£ mm 
can dominate the surface-stimulation term. 

Evaluations for the case where the nucleus of an Ih crystal is formed a) in the immersion 
mode with the basal facet on the foreign particle and b) in the contact mode with the basal 
facet on the particle and one of its prismal facet at the droplet-vapor interface can be carried 
out in a similar fashion. Besides one can consider the case where in the contact mode one 
prismal facet forms on the foreign particle and another at the droplet surface. Curiously, in 
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this situation the foreign particle does not even have to be in contact with the droplet surface. 
Moreover, the crystal nucleus may form with one of its basal facets on the foreign particle and 
the other at the droplet surface, and in this case, the foreign particle cannot be in contact with 
the droplet surface at all (unless it has a very irregular, non-compact shape). In the latter case, 
the term " contact mode" is not even appropriate. Two common features of all these "contact 
mode" situations are that: a) one of the crystal facets always forms at the droplet surface; b) 
there always exists a contact "vapor-liquid-crystal" of three water phases. Both of these factors 
(the latter even significantly stronger than the former) thermodynamically favor the formation 
of a crystal nucleus in the compact mode compared to the immersion one. The correctness 
of the term "contact mode" becomes, however, questionable, at least from a thermodynamic 
standpoint. One trivial exception from the above consideration is the case where the surface 
of the foreign particle touches the liquid-vapor interface from outside in parallel orientation 
(see Fig.l, case 4). In this situation the same facet of the crystal nucleus forms at the droplet 
surface and on the foreign particle, and there is no thermodynamic advantage for this mode 
compared to the immersion mode (when the crystal nucleus forms with the same facet on the 
same surface of the foreign particle). 

It is worth emphasizing that for accurate calculations of Wl mm and W£ on it is necessary 
to know not only the physico-chemical characteristics of the forming crystals (such as Ah, cr's 
and r's, etc..) but also the shape and size of the crystal nuclei. The latter, however, can be 
accurately determined analytically if the former are known. 

Indeed, the shape of the crystal nucleus is determined by Wulff's relations (5) and (16). For 
example, since the shape of an ice crystal cluster is known (assumed to be a hexagonal prism), 
its state is completely determined by two geometric variables (provided that its density and 
temperature are given), e.g., the height of the prism and the length of a side of a (regular) 
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hexagon (the base of the prism). However, owing to Wulff's relations, eqs.(5) and (16), only 
one of these two variables is independent. Therefore, both works W imm and W con are functions 
of only one independent variable, say, variable a, the length of a side of the hexagon. The 
concrete form of the functions W %mm = W tmm (a) and W con = W con (a) depends on the mutual 
orientation and location of the crystal cluster and foreign particle (and droplet surface in case 
of W con ). 

For instance, consider a crystal cluster formed with one of its basal facets on a foreign 
particle in the immersion mode. For the contact mode, let us consider the same basal facet 
on the foreign particle and a prismal facet (assumed to be only partially wettable by water) at 
the droplet surface. Mark the basal facets with subscripts 1 and 8 and the prismal facets with 
subscripts 2, ..., 7 (Figure 5). 

As agreed upon above, facet 8 forms on the foreign particle. Clearly, in the immersion mode 
0-O7 = a £7 = . . . = a ^, Afi = A^ = ■■■ = A^, Al 1 = A* 1 = A*" 1 . In the contact mode the 
prismal facet 7 (assumed to be only partially wettable by liquid water) represents the crystal- 
vapor interface, hence cr^ 7 = a^ 1 = ■ ■ ■ = <7g 7 , o^ 7 = <t!j . Unlike the crystal cluster in the 
immersion mode, the basal facet in the contact mode is not a regular hexagon, A^ 1 = A""' = 
A^ and A^ < A^ 7 , i.e., the surface areas of the basal facets in the contact mode is smaller 
than that in immersion mode, according to eq.(25). Let us mark two prismal facets adjacent 
to facet 7 by subscripts 2 and 6. Clearly, Af = A* 1 = A^ 1 = A* 1 , A^ = A^ = A^ < A«~< , 
Atp = Aj > Aj 1 (both inequalities are again due to eq.(25)). 

Let us use Oj and a[ (i = 2, ...,7) to denote the length of the edge formed by the basal 
facet with prismal facet i in the immersion and contact modes, respectively. In the immersion 
mode the base is a regular hexagon, i.e., a = ai = ... = 07. As clear from eq.(25), in the contact 
mode a' 2 = a' 6 < a, a' 7 > a, whereas a' 3 = a' A = a' 5 = a. 
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In the first term on the RHS's of eqs.(10) and (19) the number of molecules in the crystal 
cluster can be represented as v = p^V 1 or v = p^V 1 , respectively, where p 1 is the number 
density of molecules in phase 7 (ice) . The volume of an Ih crystal (shaped as a hexagonal prism) 
is equal to the product "height of the prism"x "surface area of the base". In both the immersion 
and contact modes the surface area of the base (regular hexagon in the former and irregular in 
the latter) is proportional to a 2 , although coefficients of proportionality are different. In both 
cases, the height of the prism, h, is linearly related to a according to Wulff 's relations (5) and 
(16), respectively. Thus, in both eqs.(10) and (19) u oc p 7 a 3 . Likewise, one can show that all 
the surface tension and line tension terms on the RHS's of eqs.(10) and (19) are proportional 
to a 2 and a, respectively. Therefore, the reversible works of formation of a crystal cluster in 
these modes can be written (tedious but simple algebra is omitted) as 

W imm (a) = -I 3 a 3 + I 2 a 2 + ha, W con {a) = -C 3 a 3 + C 2 a 2 + da, (30) 

where ^3, h, h and C3, Ci, C\ are positive coefficients, 

h = \ P ^hH®)(o-r + of - o a5 )/o«\ I 2 = ^yW + ^t ~ *"% h = 6r^ 5 , 
and 



Vs, 



a 7 7,5 _ «Sw « 7 



C 3 = ^AMn(e) {^-C h + ^ ^-{o-r + of - a a5 )W, 
C2 = « 7 + o-f ~ o a& ) (^ - C h + ^ + 

"(5 - 2(1 - (of _ a a0 )/a a 1))a ^ + {(j p _ a *p ){2 _ {o p _ ortyo**^ 1^-^ + <% - o- a «)/a«\ 
d = (t^ 5 - r Q/35 )(2 - (of - a aP )/a^) + r a "< s (5 - 2(1 - (of - a a(3 )/a^)) + 
r afiS {2C h + (2-(op -a af3 )/a^)), 

with 
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and p 1 is the number density of molecules in the crystal phase. 

Using eq.(30), one can find the length a* of a side of the hexagonal base of the crystal 
nucleus as the positive solution of the equation dW imm {a) / da\ at: = — 3I 3 a* + 2/20* + h = 0, 
or alternatively, dW con (a)/da\ at = — 3/ 3 a* + 2/ 2 o* + h = 0. which lead to a* = (2/ 2 + 



4/f + 12Ji7 3 )/6/ 3 or a* = (2C 2 + ^4C| + 12C*iC 3 )/6C 3 (see two paragraphs above eq.(30)). 
The height of the crystal nucleus (shaped as a hexagonal prism) is the same in both immersion 
and contact modes, h* = o*-^(o"^ 7 + &2 — <r )/<7p 7 . 

To numerically evaluate a* and h*, information on p, Ah,a a @), a ), a^ 7 , cjp 7 ,a" 7 , cr^ 7 ) p , 
r "P7 ; r °P" ^ r «7o ; anc j r P7 5 i s needed. Experimental data on p 7 , A/i, and a a ^) are readily available 
(even as functions of temperature). For our evaluations they were taken to be p" f = 0.92^^/18 
cm~ 3 , Ah ~ 333.55 x 10 7 Na/18 erg (where Na is the Avogadro constant), and a a @) = 83 
dyn/cm. Some data on a^ 1 and u" 7 as well as on the mean value of the crystal-vapor surface 
tension have been also reported (see ref.2 for a short review). However, virtually no reliable 
data are currently available for any solid-ice interfacial tensions and line tensions in "solid 
substrate-ice-liquid water-water vapor" systems. These were hence chosen somewhat arbitrarily, 
the main criterion being a reasonable agreement of the estimates extracted from the above 
equations with those obtained from the experimental data on homogeneous ice nucleation. 
Considering nucleation of ice crystals on a foreign particle such that at a given temperature 
(j a ° = 40 dyn/cm, cr^ 7 = 23 dyn/cm, a" 7 = 24 dyn/cm, a 7 = 50 dyn/cm, & )p = 
102 dyn/cm, T a ^ = -10~ 4 dyn, T af3& = 7 x 10" 5 dyn, r" 7 * = 10~ 5 dyn, and t^ 5 = 
5 x 10~ 6 dyn (reasonable choice according to scarce data available in literature), equations for a* 
and h* would provide a* ~ 29x 10~ 8 cm, h* ~ 35x 10" 8 cm, and W™ n -W\ mm ~ -12k B T ht . As 
intended, these values are consistent with the estimates obtained above from the experimental 
data on the homogeneous nucleation rate. 
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4 Concluding Remarks 

Previously, in the framework of CNT a criterion was found for when the surface of a droplet can 
stimulate crystal nucleation therein so that the formation of a crystal nucleus with one of its 
facets at the droplet surface is thermodynamically favored (i.e., occurs in a surface stimulated 
mode) over its formation with all the facets within the liquid phase (i.e., in a volume-based 
mode). For both unary 21 and multicomponent 22 a droplets, this criterion coincides with the 
condition of partial wettability of at least one of the crystal facets by the melt (the contact 
angle, measured inside the liquid phase, must be greater than zero). 

However complex a theory of homogeneous crystal nucleation in droplets may be, the pres- 
ence of foreign particles, serving as nucleating centers, makes the crystal nucleation phenomenon 
(and hence its theory) even more involved. Numerous aspects of heterogeneous crystal nucle- 
ation still remain obscure. One of most intriguing problems in this field remains the strong 
enhancement of heterogeneous crystallization in the contact mode compared to the immersion 
one. It has been observed that the same nucleating center initiates the crystallization of a 
supercooled droplet at a higher temperature in the contact mode (with the foreign particle just 
in contact with the droplet surface) compared to the immersion mode (particle immersed in the 
droplet). 2 Many heterogeneous centers have different nucleation thresholds when they act in 
contact or immersion modes, indicating that the mechanisms may be actually different for the 
different modes. Underlying physical reasons for this enhancement have remained largely un- 
clear, but the phenomenon of surface-stimulated (homogeneous) crystal nucleation had strongly 
suggested that the droplet surface could enhance heterogeneous nucleation in a way similar to 
the enhancement of the homogeneous process. 

In this paper we have extended the approach, previously developed in refs.21,22, to hetero- 
geneous crystal nucleation on a solid particle (in both immersion and contact modes) and have 

25 



presented a thermodynamic model shedding some light on the mechanism of the enhancement of 
this process in the contact mode. Our thermodynamic analysis suggests that the droplet surface 
can indeed thermodynamically enhance crystal nucleation in the contact mode compared to the 
immersion mode. Whether this occurs or not for a particular foreign particle is determined, 
however, by the interplay between various surface tensions and four line tensions involved in 
this process. As clear from our model, the droplet surface may stimulate the heterogeneous 
crystal nucleation even in the case where the foreign particle is actually completely immersed 
therein, but is situated closely enough to the surface. This suggests that the term "contact 
mode enhancement" is probably not very appropriate for this phenomenon. 

As a numerical illustration of the proposed model, we have considered heterogeneous nu- 
cleation of Ih crystals on generic macroscopic foreign particles in water droplets at T = 253 
K. Our results suggest that while the droplet surface always stimulates crystal nucleation on 
foreign particles in the "contact mode" , the line tension contribution to this phenomenon (due 
to the contact of three water phases, "vapor- liquid-crystal" ) may be as important as the surface 
tension contribution. 
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Captions 

to Figures 1 to 5 of the manuscript 

"Thermodynamics of heterogeneous crystal nucleation in contact and immersion 

modes" by Y. S. Djikaev and E. Ruckenstein 

Figure 1. Heterogeneous crystal nucleation in a liquid droplet surrounded by vapor. "Im- 
mersion" mode: the crystal cluster forms with one of its facets on a foreign particle, completely 
immersed in a liquid droplet;ia all other crystal facets interface the liquid. "Contact" mode: the 
foreign particle is in contact with the droplet surface; the cluster forms with one of the crystal 
facets on the particle, another facet at the liquid-vapor interface, and all other facets making 
the "crystal-liquid" interface. Cases 1 through 4 represent a few of possible variations of the 
"foreign particle-droplet surface" contact. 

Figure 2. Illustration to Wulff's relations. The surface area and surface tension of facet i are 
denoted by Ai and <7j, respectively; hi is the distance from facet i to reference point O. 

Figure 3. Heterogeneous formation of a crystal nucleus on a foreign particle completely im- 
mersed in the liquid. 

Figure 4. Heterogeneous formation of a crystal nucleus on a foreign particle in contact with the 
liquid- vapor interface. 

Figure 5. Heterogeneous formation of an Ih cluster on a foreign particle in the immersion and 
contact modes, contact with the liquid-vapor interface. The crystal cluster has a shape of a 
hexagonal prism. One of the basal facets (facet 8) is formed on the foreign particle, the other 
(facet 1) interfaces the liquid. Two of the prismal facets (with numbers 4 and 7) lie in the plane 
perpendicular to the Figure. Prismal facets 5 and 6 cannot be seen by the reader, so they are 
shown in the parentheses. 



30 




Vapor, phase p 



Figure 1: 



31 



Vapor 



Liquid 




Figure 2: 



32 



Foreign particle 

(phase 8) 



x x 



Liquid 

(phase a) 




Crystal 
(phase y) _ ' 



/A o<*1 



Figure 3: 



33 



Vapor 
(phase P) 

X.-1 ° 




Foreign / 

particle / \ 
(phase 8) / V s *' 

/ Liquid 


\ 


■' O 

Crystal 
(phase y) 


^\y (phase a) 







Liquid 

(phase a) 



Figure 4: 



34 



Vapor 

(phase p) 




Vapor 

(phase p) 



Liquid 

(phase a) 




Immersion mode 



Liquid 
(phase a) 



2(6) 
3(5) 



Contact 



Liquid 
(phase a) 



mode 



Crystal 

(phase y) 



Figure 5: 



35 



